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ABSTRACT 

The star XO-5 (GSC 02959-00729, V=12.1, G8V) hosts a Jupiter-sized, i? p =1.15±0.12 R, 7 , tran- 
siting extrasolar planet, XO-5b, with an orbital period of 4.187732±0.00002 days. The planet mass 
(M p =1.15±0.08 Mj) and surface gravity (g p =22±5 m s -2 ) are significantly larger than expected by 
empirical M p -P and Af p -P-[Fe/H] relationships. However, the deviation from the M p -P relationship 
for XO-5b is not large enough to suggest a distinct type of planet as is suggested for GJ 436b, HAT- 
P-2b, and XO-3b. By coincidence XO-5 overlies the extreme H I plume that emanates from the 
interacting galaxy pair NGC 2444/NGC 2445 (Arp 143). 

Subject headings: planetary systems - stars: individual (GSC 02959-00729) - galaxies: individual 
(NGC 2444, NGC 2445) 
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1. INTRODUCTION 

We report the discovery of a Jupiter-sized planet, XO- 
5b, that transits the G8V, V=12.1, star GSC 02959- 
00729 (XO-5) with an orbital period, P~ 4 days. 
XO-5b provides a valuable addition to the empirical 
pop ulation trends amo ngst transiting extrasolar plan- 
ets (jTorres et al.ll2008T ). especially at P>4 day periods 
where only four other well characterized transiting plan- 
ets are known (O GLE-TR-lllb, HAT-P-lb, HAT-P-2 b, 
fc HP 17156b. iPont et all l200i iBakos et all [20030; 
iBarbieri et al.ll2007l respectively). With a sampl e of the 
first six known transiting planets, iMazeh et al.l (|2005D 
pointed out the possibility of a linearly decreasing trend 
of pla net mass as a func tion of increasing orbit al pe- 
riod. iTorres et all (|2008l ) and ISouthworthi (|2OO80 have 
recently revisited this and other relationships amongst 
the now larger sample of transiting planets analyzed in a 
homogeneous fashion. As shown in Section 2J XO-5b is 
one of the few planets that is inconsistent with general 
trends based on the currently known sample of transiting 
Jupiter-sized planets. 

In Section [2] we describe the photometric and spec- 
troscopic observations that we used to discover and 
characterize the XO-5, XO-5b system. Section 13.11 de- 
scribes our spectroscopic determination of the atmo- 
spheric properties of XO-5, which we combine with the- 
oretical isochrones to constrain the stellar mass and ra- 
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dius. The mass estimate of XO-5 is employed along with 
a high precision transit light curve to measure the radius 
of XO-5b in Section 13.21 Radial velocity measurements 
determine the mass of XO-5b as outlined in Section l3~3l 
In Section l3~4l the database of transit light curves enables 
refining the transit ephemeris for XO-5b. We conclude 
in Section [4] how the discovery of XO-5b adds to the un- 
derstanding of the population of transiting Jupiter-sized 
planets. 

2. OBSERVATIONS 
2.1. XO Survey Photometry 

iMcCullough et all (|2005l ) provides details of the 
XO survey; additional details of the candidate tran- 
siting planet select io n pr ocess can be found in 
IMcCullough fc Burkd ([27)01 . XO-5b is the fifth 
Hot Jup iter transiting pla n et an n ounced by t h e XO 
survey (IMcCullough et all [20061 iBurke etjap [20071 : 
I Johns-Krull eTafl [20081 : IMcCullough eTaTl l2T)f)l . and 
was identified as a candidate from the 62° x 7° field of 
view centered on 8 h RA. This field was observed at a 10 
min cadence with an overall rms noise, a — 0.011 mag, 
for this V=12.1 star during two seasons: Nov. 2003-Mar. 
2004 and Nov. 2004-Mar. 2005. XO-5 is well isolated; 
XO-5 contributes 90% of the flux in the 75" radius aper- 
ture used for XO photometry. 

Figure [T] shows the XO light curve phased to the pe- 
riod with the most significant transit even t identified by 
the B ox-fitting Least Squares algorithm ([Kovacs et alJ 
I2002D . In the phased light curve, two full transit events 
and five partial events contribute to the detection. After 
subtracting a transit model (determined from the high 
precision light curve in Figure [3]), the average residual 
RMS during transit nights is a = 0.008 mag (30% lower 
noise than the average for the entire light curve). Adding 
in quadrature the Signal to Noise Ratio (SNR) for each 
transit yields an overall detection SNRrotaZ ~ 12. It is 
not routine to follow up all candidates with this low of 
SNR^ota;- However, SNRr ot q; is only one of several crite- 
ria for selecting ca ndidates ([McCullough fc Burkd 120071 : 
IBurke et al.l 120071 ). The efficiency of photometric fol- 
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Fig. 1. — Discovery phased light curve from the XO Project data. 
(Top) The phased light curve around the transit event. (Bottom) 
The phased light curve over the full orbital period. The individ- 
ual measurements (points) are shown binned (solid line) to reduce 
noise. The transit occurs at At = 0.0 day and any secondary 
eclipse for a circular orbit would occur at At = 2.1 day (dotted 
line). No secondary eclipse is evident above the noise. 

lowup provided by the XO Extended Team is an addi- 
tional resource of the XO project to aid in following up 
these lower SNR/r ta/ candidates. 

2.2. XO Extended Team and Follow Up Photometry 

The Extended Team (E.T.) provides photometric fol- 
low up of XO candidates. The E.T. (M. F., C. F., 
E. G-M., J. G., F. M., G. M., and T. V.) is a 
collaboration of professional a nd amateur astronomers 
(|McCullough et all 120051 120061 ) . The E.T. obtains light 
curves of XO candidates at higher angular resolution and 
guides the photometric and spectroscopic follow up nec- 
essary to classify a candidate as a bona fide planetary 
companion. Table Q] provides E.T. photometry for XO-5. 
For the E.T. light curves, the median differential magni- 
tude out of transit provides the flux normalization and 
the standard deviation out of transit provides the un- 
certainty in the measurements. Figure [2] shows the E.T. 
light curves obtained for XO-5b. 

On 2008 January 20 we observed a transit event of 
XO-5b with the 1.8-m Perkins Telescope at Lowell Ob- 
servatory using th e PRISM instrument in imaging mode 
(| Janes et al.ll2004l ). The observations arc a continuous 
series of R-band images covering a 10.8' x 8.8' subar- 
ray field of view with 10s exposure time and 30s readout 
resulting in an overall t ca d =40 s cadence. The detec- 
tor has 0.36" pixels, and the FWHM seeing was 1.8". 
The target was observed over the airmass range, 1.12 < 
X < 1.6, with observations commencing at high airmass. 
The evening was nonphotometric; several percent trans- 
parency variations in raw photometric magnitudes ex- 
ist around the linear trend with airmass. To remove 
the transparency variations a differential light curve is 
calculated employing four other comparable bright stars 
as comparison. The resulting differential light curve of 
XO-5 is shown in Figure [3] The out of transit data are 
limited in duration. Thus, minimizing the scatter dur- 
ing the flat, in-transit portion of the light curve, guided 
the choice of photometric aperture size and comparison 
stars. The resulting noise during the in-transit portion 
is a. 
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0.002 mag for each exposure. This is ~ 20% 
higher noise than e xpected from Poisson, scintillation 
(|Dravins et al.|[l997t ). sky, and read noise. More relevant 
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Fig. 2. — Binned light curves from the Extended Team for XO- 
5. The text in color matches the corresponding light curve and 
indicates the date and passband of the observations. 
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Fig. 3. — Top: Light curve from the 1.8m Perkins Telescope at 
Lowell Observatory in the R-band (points), along with the best-fit 
transit model in a x 2 sense during the MCMC analysis (solid line). 
Bottom: Residual of data from the best-fit transit model. 

for transit detection and characterization is the photo- 
metric noise rate, a pnr = o- exp yJt ca( i/QQs = 0.0016 mag 
mn~ | which takes into account the dead time of detector 
readout. The calibration and photometry was performed 
in IRAF with a customization in the IRAF phot task to 
report magnitudes to 10 -4 precision, rather than the de- 
fault of 10~ 3 precision. 

We obtained photometric B, V, Rc, and Ic magnitudes 
for XO-5 using a 0.35-m telescope (Tabled]). Photomet- 
ric observations were obtain ed on three se parate nights 
with 2-6 Landolt star fields (|Lan dolt 1992|) interspersed 
with the XO-5 field in order to define the color and air- 
mass transformation from the instrumental system to the 
standard system. Each night 20-40 Landolt stars defined 
the transformation with rms scatter of 0.04, 0.03, 0.02, 
0.03 for the BVRI passbands, respectively. The result- 
ing standard photometry in Table [5] is the weighted aver- 
age amongst the thre e nights. Standard star photometry 
of XO-5 from TASS (jDroege et a l.||2006|) agrees with our 
photo metry. Photometry from 2MASS (jSkrutskie et al.1 
I2006D is also given in Tablelg] XO -5 is saturated in SDSS 
( Adelman-McCarthv et al.l2008D images and too faint for 
the TYCHO-2 (\R&e et al.ll2000li database. 

2.3. Arp 143 - Potential Archival Measurements 

By coincidence, XO-5 lies near the line of sight to the 
interacting galaxy pair NGC 2444/5 cataloged as Arp 
143 in the 'Material Emanating from E Galaxies' group 
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Fig. 4. — False color image for finding XO-5. The box encloses 
XO-5, and the scale bar indicates 1 arcmin. The finder is ~9 
arcmin per side, and North is up and East is toward the left. The 
interacting galaxy pair Arp 143 is visible in the lower left. 
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Fig. 5. — Digital Sky Survey image (Red band) of XO-5 (arrow) 
showing the coincidental H I plume (contour) that emanates from 
the interacting galaxies Arp 143. The H I em ission contours are 
VLA observations from Applcton et al. (1983). 



(|Arp| [l966'). Figure [4] shows a false-color finder chart for 
XO-5 enclosed by the box along with the Im peculiar 
galaxy NGC 2445 (patchy blue galaxy) and smoother SO 
peculiar galaxy NGC 2444 three arcmin South of XO-5. 
The VLA has detected one of the largest known (14.5' 
in extent) H I plumes, which emanates from Arp 143 
that directly coincides with XO-5 as shown in Figure EJ 
(lAppleton et al.lll987l : Iffigdon et allll997t iHibbard et all 
120011 ). With a velocity shift of 4000 kms" 1 , it is improb- 
able that the H I emission from Arp 143 can be attributed 
to XO-5. Archival observations of Arp 143 may provide 
long term photometric monitoring or serendipitous tran- 
sit data for XO-5, possibly constraining ephemcris vari- 
ability. 

2.4. Spectroscopy 

After confirmation of the XO transit light curve from 
E.T. observations (see §2.2|) . we initiated queue sched- 
ule observations of XO-5 with the High- Resolution Spec- 
trograph (HRS), a fiber fed cross-dispersed cchcllc spec- 



trograph (jTullI 119981) . on the 11-m Hobby- Eberly Tele- 
scope (HET) located at McDonald Observatory in or- 
der to measure the mass of the planet. HRS observa- 
tions with an iodine gas cell to yield precise radial ve- 
locities commenced on 2007 December 7. Table [3] pro- 
vides dates of the HRS observations along with the re- 
sulting radial velocities. The instrument setup provides 
R=63,000 resolution and covered the wavelength range 
4000 < A < 7800 A, centered at A = 5900 A. The two- 
dimensional echel le spectra are extra cted using proce- 
dures described in lHinkle et all (|2000f ). Radial velocities 
for XO-5 are determined in § 13.31 

To measure the stellar parameters of XO-5, we also 
obtained HRS HET spectra without the iodine absorp- 
tion cell, using the 5936 A and 6948 A settings of the 
cross-disperser to cover the wavelength intervals 5150 
< A < 5200A and 6000 < A < 6200A at R=63,000. 
Three spectra at each configuration were co-added for 
analysis. The SNR is 35 and 50 for the co-added spec- 
trum in the blue and red configurations, respectively. 
The XO-5 stellar parameters are determined in § 13.11 

3. ANALYSIS 

3.1. Stellar Properties 

A transit light curve and radial velocity measure- 
ments do not completely solve the system of stellar and 
planet propert ies; these measurements deter mine the 
stellar density (jSeager fc Mallen-Ornelasll2003l ). Solving 
for the stellar and planet properties independently re- 
qui res an additional c onstraint. Following the procedure 
of (jBurke et al.ll2007l ) in analyzing the transiting planet 
XO-2b, spectroscopic analysis combined with stellar 
isochrones provides an estimate of the stellar mass, M*, 
and its uncertainty. This prior constraint on M 4 enables 
solving the system completely. We u sed the Spectroscopy 
Made Easy (SME) analysis package (IValenti fc Piskunovl 
119961 ) with refinements from IValenti fc Fischer! (|2005f ) to 
determine the spectroscopic the properties of XO-5. 

Table lists XO-5 stellar parameters from the SME 
analysis of a single coadded spectrum. The 1-a uncer- 
tainties in the stellar parameters given in Table [2] are 
based on the typical rms scatter in parameters mea- 
sured in independent, multiple spectra for stars in the 
SPOCS catalog. XO-5 has an enhanced metal abun- 
dance, [Fe/H]=0.25±0.03, and the effective tempera- 
ture, T e ff and surface gravity, logg are consistent with 
a G8V spectral type according to Appendix B of iGravl 
(|1992l ). The empirical Te.ff, intrinsic color calibra- 
tions of IWorthev fc Led (|2006f ) predict a color excess, 
E(B - V) = 0.03 ± 0.04 and E(V - K) = 0.03 ± 0.03 for 
XO-5. This difference is negligible within the uncertain- 
ties, and a typical reddening law predicts E(V —K) prec i = 
2.74 x E(E> — V) = 0.08, lower tha n the measured excess 
E(V - K). The extinction map of ISchlegel et all (|1998f ) 
also measures a modest E(B — V) = 0.05 toward XO-5. 
During the analysis we assume that the standard pho- 
tometry is free of reddening. 

Using the primary observables from the SME analy- 
sis {T e ff, abundances, and logg) and the apparent V- 
band magnitud e (§I2.2 |), we de termined M* and i?* from 
Y 2 isochrones (lYi et all 120011 ). using the procedure of 
IValenti fc Fischer! (|2005f ) . M+ and R+ depend on the un- 
known distance to XO-5, thus the probability density 
function for M*, i?*, and age are calculated for a se- 
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Fig. 6. — Parametric curves for the Radius (Top) and Mass 
(Bottom) of XO-5 as a function of the surface gravity from the SME 
isochrone analysis, respectively (see § 13. It . The SME spectroscopic 
determination of surface gravity for XO-5 (solid), the l-cr internal 
uncertainty (dashed), and the l-cr systematic uncertainty (dotted) 
are shown with vertical lines. 

quencc of trial distances in steps of 10 pc. Estimates 
of M* and i?+ from this isochrone analysis provide an 
additional estimate of log<7i SO as a function of distance 
to XO-5. Figure [5] shows R+ and M* as a function of 
logfeo- A lower mass, lower luminosity, younger age, 
and smaller distance star for XO-5 corresponds to higher 
surface gravity (left side of Figure 0; whereas, a higher 
mass, higher luminosity, older age, and larger distance 
for XO-5 corresponds to lower surface gravity (right side 
of Figure [6]). 

A primary goal of the SME analysis is to provide an 
estimate of M* and its uncertainty in order to provide 
a unique solution when fitting the light curve and ra- 
dial velocity data. The spectrum alone yields an esti- 
mate of the stellar gravity, logc/ sme that is independent 
of \oggi S0 . The condition logg sme =loggj SO and the un- 
certainty in log<7 sme directly yields M* from the data 
shown in Figure [6l The solid vertical line in Figure [6] 
shows \ogg sme which corresponds to M* = 1.0M©. The 
long dashed vertical line in Figure [6] represents the uncer- 
tainty of logg sme , aig CJ = 0.06 dex. This uncertainty is an 
internal precision of objects in the SPOCS catalog ana- 
lyzed in a uniform manner. For the purpose of estimating 
Mi,, a more appropriate uncertainty in log<? is determined 
through comparison of logg sme from the SPOCS catalog 
to independent logq me asurements from the literature. 
IValenti fc Fischer! (|2005[ ) find rms scatter o~i gg = 0.15 dex 
by comparing the SPOCS catalog to external catalogs. 
The dotted vertical line in Figure [6] shows this larger un- 
certainty in logg. M* varies weakly with logg, and we 
adopt an uncertainty of, er M = O.O3M for M*, but also 
explore the impact a larger, cxm = 0.07 Mq, has on the 
system parameters. 

3.2. Markov Chain Monte Carlo Light Curve Analysis 

iFordl (|2005l ). iGregorvl (|2005f ). and references therein 
provide a thorough discussion of the theory behind 
Markov Chain Monte Carlo (MCMC) Bayesian analysis 
along with a practical MCM C implementation for radial 
velocity planet detection, and lHolman et all (|2006f ) intro- 
duced MCMC techniques to the analysis of a transiting 
planet light curve. To analyze the X O-5b light curve, we 
follow the MCMC implementation of iBurke et~aH (|2007f ) 
where they analyzed the light curve of the transiting 



planet XO-2b. 

To calculate the Bayesian posterior probability for the 
system parameters, the likelihood function is given by 
e -o.5x ) w here we have assumed the errors are normally 
distributed, and the data have uniform weights, x 2 is the 
squared differenc e between observations and the analytic 
transit model of iMandel fc Agoll (|2002l ). The model as- 
sumes negligible eccentricity. The observations are the 
R-band data from the 1.8m Perkins Telescope shown in 
Figure [3l To provide the additional constraint necessary 
to uniquely determine the planet properties, we adopt 
an informative prior on M+ that is a Gaussian with cen- 
ter M+ = I.OMq and standard deviation a = O.O3A/0 
(see § I3.ip . We adopt uninformativc priors for R+, R p , 
and i that are uniform. To improve the efficiency of the 
MCMC calculation, the computation is done using the 
set of parameters i?*, p — R p /R±, and the total transit 
durati on from 1 st to 4 th contact, r instea d of R+, R v , 
and i ([Burke et al.ll20"07h . Appendix A of IBurke et all 
(|2007| ) provides the form of the priors necessary to main- 
tain the uniform priors in R±, R p , and i when the cal- 
culation is done using i?*, p, and r. The limb darkening 
coefficients, u\ and U2, of the transit model are allowed 
to vary. Limb darkening is described by a quadratic law, 
I = 1 — Wi(l — /i) — 112(1 — p) 2 , where / is the specific 
intensity normalized to unity at the center of the stellar 
disk and p is the cosine of the angle between the line 
of sight and t h e surf ace normal. In practice, we follow 
iHolman et al.l (|2006f ) by adopting a\ = u\ + 2u 2 and 
<22 = 2ui — U2 as the parameters used in the calcula- 
tion. The limb darkening coefficients are constrained by 
requiring the highest surface brightness to be located at 
the disk center (u\ > 0.0), by requiring the specific in- 
tensity to remain above zero (u\ + 1*2 < 1.0), and by not 
allowing limb-brightened profiles (ui+2m2 > 0.0). The fi- 
nal free parameter is the mid-transit time offset, t , from 
the initial ephemeris with a period obtained from ana- 
lyzing the XO observations and E.T. observations (see 

The Markov Chain analysis employs the Metropolis- 
Hastings algorithm with a normal proposal distribution 
and a multiple block sampling technique where each step 
in the chain consists of a number of intra steps updating 
each individual parameter in turn. Several short, trial 
chains iteratively yield scale factors of the normal pro- 
posal distribution for each parameter with a 25% to 40% 
acceptance rate for the trial samples. In addition to the 
choice of i?*, p, and r as the variables for the calcula- 
tion, a linear transformation between parameters yields 
an eigenbasis set of parameters with a multi-normal non- 
covariant relationship that further r educes the mutua l 
degeneracy amongst the parameters (|Burke et al.ll2007t) . 
These transformations yield an efficient MCMC calcu- 
lation with an autocorrelation length N cor = 5 for the 
chain. The estimate of the posterior probability comes 
from 7 independent chains of length N c h n = 70000 with 
varying initial conditions. This results in an effective 
length Net, = N chn /N cor = 14000. 

Figure \J\ shows the resulting posterior probability dis- 
tribution for select system parameters after marginaliza- 
tion over the other parameters. The posterior proba- 
bility is simply a normalized histogram of the MCMC 
sample values. We adopt the median as the best sin- 



gle point estimate of the posterior probability. The \-a 
credible interval for a parameter is given by the symmet- 
rical interval around the median that contains 68% of 
the samples. Figure [8] shows the joint posterior proba- 
bility for the stellar limb darkening coefficients, u\ and 
U2- The constraints on the limb darkening coefficients are 
not very constraining. Figurc[S]shows theoretically calcu- 
lated limb dark ening coeffici ents in various photometric 
passbands from lClaretl (2000) for a star with the physical 
properties of XO-5. The best estimates for limb darken- 
ing coefficients (diamond symbol) from the R-band tran- 
sit light curve is consistent with the theoretical R-band 
coefficient (triangle symbol). 

Using the SME isochrone analysis, the MCMC sam- 
ples for R± translate into distance to and age estimates 
of XO-5. Similar to the procedure that defines the prior 
on M*, The SME isochrone analysis provides a para- 
metric relationship between i?* versus trial distance and 
age (i.e., for a given stellar radius estimate, the SME 
isochrone analysis has a best distance and age estimate 
for the system). 

The theoretical planet models of iFortnev et al.l (j2007l ) 
characterize the amount of heavy elements (elements 
heavier than H & He) for a given orbital semi-major axis, 
R p , and M p . The best estimate of R p for XO-5b is consis- 
tent with zero heavy element content, and the 1-a upper 
limit is 5 Mq of heavy elements. Table [5] summarizes the 
properties of XO-5b. 

The uncertainties given in Table [5] for XO-5b and Ta- 
ble [2] for XO-5 represent the internal precision of the 
experiment. These uncertainties represent the expected 
scatter of values obtained if the experiment was repeated 
with similar quality data and identical procedures. Other 
systematic sources of error are most likely comparable or 
larger than the internal precision. The sources of sys- 
tematic error only enter into the prior for M*. In light 
of these potential sources of systematic uncertainty, the 
analysis is repeated with an increased standard deviation 
of the Gaussian prior on M* to a = 0.07 M Q . This larger 
uncertainty on M* did not increase the uncertainties in 
the other parameters beyond the most significant digit 
as given m Tableland Tabled This indicates the pho- 
tometric noise rather than uncertainty in dominates 
the uncertainty in the system parameters. 

After fixing the prior on Af*, the light curve analy- 
sis yields an additional estimate of \ogg rncmc as given in 
Table El logg mcmc does not differ significantly from the 
SME based logg sme . We redetermined the stellar prop- 
erties by fixing \ogg sme =\ogg mcmc when analyzing the 
spectra. This reduced M* = 0.02 M Q . This variation 
is within our original uncertainty of M*. Thus, we did 
not iterate the MCM C analysis with the SME analysis 
(jSozzetti et al.ll2007t) . 

3.3. Radial Velocity Measurements 

We determined the mass of XO-5b using the radial ve - 
locity techniques describ ed bvlMcCullough et all (|2006D , 
iBurke et al.1 (|2007t ). and Uohns-Krull et al.l (|2008h . We 
summarize the procedure here, where we derive radial 
velocity information from the HET spectra (see § 12. 4p . 
An iodine gas cell imprints iodine absorption lines on 
the spectrum of XO-5 and this is compared to a model 
spectrum to obtain radial velocity shifts with respect to 
the topoccntric frame. We construct model spectra by 




Fig. 7. — Marginalized posterior probability for the XO-5 and 
XO-5b parameters from MCMC samples. We adopt the median 
of the posterior probability (solid) vertical line as the point esti- 
mate of the parameter. The dashed lines indicate the 68% credible 
interval. 




Fig. 8. — Joint posterior probability for the stellar limb darken- 
ing coefficients, u\ and ui. The solid contours are isoprobability 
contours containing 68% and 90% of the MCMC samples. The 
remaining 10% of the samples lying outside the region of highest 
probability are also shown (points). The theoretically calculated 
limb darkening coefficients are shown for various passbands (trian- 
gle) along with the best estimate limb darkening coefficients from 
the MCMC analysis of the R-band transit light curve (diamond) . 
The R-band light curve coefficients are consistent with the theoret- 
ically calculated R-band coefficients. The prior limits for the limb 
darkening coefficients are indicated with dashed lines. 

multiplying a very high-resolution FTS spectrum of the 
Sun (scaled to better match the observed line depths of 
XO-5) times a very high-resolu tion FTS spectrum of an 
iodine gas cell ([Coch ran 2 000T) and then convolving the 
result with a Voigt profile to approximate the line-spread 
function of the instrument. A radial velocity shift is de- 
termined independently for each ^15 A row of the 2-D 
echelle spectrum of XO-5 over the wavelength range with 
strong iodine absorption, 5210 < A < 5700 A. The indi- 
vidual stellar radial velocity estimates from each of the 
15 A sections (after transformation to the barycentric 
frame) determines the stellar radial velocity measure- 
ment at each epoch and its associated 1-a uncertainty 
(see Table [3]). 

Figure [9] shows the resulting radial velocity curve 
phased with the XO-5b ephemeris determined from the 
transits and assuming zero eccentricity. The typical un- 
certainty per epoch is a^v = 20 — 40 m s _1 , limited 
mainly by Poisson noise in the spectrum. The radial 
velocity semi-amplitude, K=145±10 m s . This ampli- 
tude results in M p =I.15±0.08 Mj for XO-5b, assuming 
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Fig. 9. — a) The radial velocity of XO-5 oscillates sinusoidally 
with a semi-amplitude K = 145±10 m s _1 , implying XO-5b's mass 
is 1.15±0.08 Mj. b) The period and phase of the radial velocities 
were fixed at values determined by the transits. The mean stellar 
radial velocity with respect to the solar system's barycenter has 
been subtracted. In order to determine K, we used the HET spectra 
calibrated with an iodine absorption cell (filled circles), c) In this 
representation of the data, a circular orbit yields a straight line of 
slope — K. 

M*=1.0 M Q for XO-5 and a circular orbit for XO-5b. 

3.4. Ephemeris and Transit Timing Variations 

The XO, E.T., and Perkins Telescope transit light 
curves enable refining the ephemeris for XO-5b. The 
high precision light curve from the Perkins Telescope 
(see § I2.2[) . provides the estimate for the mid-transit zero 
point of the ephemeris as well as the transit model cm- 
ployed to derive mid-transit timing for the other lower 
precision light curves from the E.T. and XO survey. 
The MCMC samples from the Perkins Telescope light 
curve analysis (lower right panel of Figure [7| provide 
the best estimate and error for the ephemeris zero point, 
HJD =2454485.6664±0.0004 (a = 35 s). 

To refine the orbital period of XO-5b, the remaining 
XO and E.T. light curves (E.T. light curves are shown in 
Figure^) provide mid-transit timing estimates. The best 
fitting transit model in a \ 2 sense to the high precision 
Perkins Telescope light curve (Figure [3]) is fit to the other 
light curves with only mid-transit time and out of transit 
flux zero point as the only free parameters. The expected 
variations in the transit model as a function of photomet- 
ric bandpass is negligible compared to the photometric 
noise in the E.T. and XO light curves, so the limb dark- 
ening coefficients are fixed to the values from analyzing 
the R-band Perkins light curve. Table [4] provides mid- 
transit timings along with their uncertainty. The XO 
observations have 10 min cadence and ~1% photometric 
noise making transit timing uncertain to 28 min. This 
uncertainty comes from comparing the a transit model fit 
a by-eye estimate of the mid-transit timing. Within the 
uncertainties there is not a significant systematic differ- 
ence between mid-transit timings estimated by-eye and 
from the model fit. The much higher quality E.T. light 
curves have an uncertainty of 4.3 min in mid-transit tim- 
ing. This uncertainty is estimated from comparing the 
groups of mid-transit timing observations for the same 
event but different observers. 

With these uncertainties, an ephemeris model with pe- 
riod as the only free parameter results in \ 2 = 12-6 with 
v = 14 degrees of freedom indicating the timings are con- 
sistent with a fixed period, P =4.187732±0.00002 day. 



The ephemeris for XO-5b accumulates a 5 min uncer- 
tainty by 2010. Given the large uncertainty in transit 
mid-point timing from the XO survey data, the period 
was calculated using E.T. data alone. The earliest, and 
only, E.T. light curve from the 2006-2007 observing sea- 
son for XO-5b provides the strongest constraint on the 
orbital period. The remaining E.T. light curves are from 
the 2007-2008 observing season for XO-5b. The period 
derived from E.T. data alone is within l-u of the deter- 
mination that includes the XO survey data. 

4. DISCUSSION & CONCLUSION 

The transit candidate of Mandus hev et al] (|2005l ) il- 
lustrates the non-negligible potential for triple stars to 
have transit light curves and radial velocity variations 
consistent with a planet. However, in the case of XO-5b, 
attempts to explain the light curve and spectroscopy with 
a physical stellar triple fail. We employ the Y 2 isochrone 
appropriate for the the physical properties of XO-5 sup- 
plemented with the low- mass stellar isochrone between 
0.072< M* <0.5 M from lChabrier et all (120001) . stellar 



limb darkening coefficien ts fromlClaretl ( 2000D . the light 
curve synthesis routine of lWilsonl (|1993f ) , and assume any 
potential stellar binary has zero orbital eccentricity to 
model a stellar triple system. The constraints on the 
transit duration and transit depth from the light curve 
require M* > 0.90 M Q for the primary of a stellar binary 
blended with the light of XO-5. The required primary 
of the stellar binary has >45% the flux of XO-5 and 
has a radial velocity semi-amplitude, K > 18 km s . 
Such a binary would be readily apparent in the spec- 
trum of XO-5 given the narrow spectral features for XO- 
5, wsini < 3.3 km s _1 . We cannot completely rule out 
the possibility of a line-of-sight faint background binary 
blended with the light of XO-5 as an explanation for the 
observations. However, the sinusoidal shape of the radial 
velocity variations necessitate the line-of-sight binary to 
have a systemic velocity similar to XO-5 otherwise the 
radial vel ocity curve develop s asymmetries that are not 

o bserved ([Torres et al.ll2005l). 

iTorres et al.l (|2008t) and iSouthworthl (|2008f ) recently 



derived transiting planet properties in a uniform fash- 
ion in order to further test various trends amongst the 
transiting planets that have previously been examined 
with smaller and less homogen eous samples of transit- 
ing planets (jMazeh et al.l I2005D . T he la rger sample of 
transiting planets in ITorres et al.l (|2008l ) still shows a 
general trend for decreasing M p with increasing orbital 
period, but there is significant scatter and some very sig- 
nificant o utliers (GJ 43 6 b, HA T-P-2b, and XO-3b). Fur- 
thermore, ITorres et all (|2008l ) find a M p -P-[Fe/H] rela- 
tion that reduces the scatter in the M p -P relation. At 
fixed period, the more metal poor stars host higher mass 
planets. XO-5b has AM p = 0.6± 0.08 M, T high e r mas s 
than then M p -P trend defined in ITorres et all (|2008D . 
The metallic i ty cor rection to the M p -P trend found by 
ITorres et all (|2008t) . AM p = -0.14 Mj in the case of 
XO-5b, is o pposite sign f r om w hat is measured. As com- 
mented bv iTorres et all (|2008fl , the M p -P-[Fe/H] rela- 
tion does have a scatter larger than observational uncer- 
tainties. Thus, additional physical conditions impact the 
M p -P relation, the full variety of Jupiter-class extrasolar 
planets has not been fully explored, or the trend is more 
relevant to shorter orbital periods than XO-5b. For in- 
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stance, iMazeh et all (|2005[) propose the M p -P relation is 
setup by the thermal evaporation of planets too small to 
survive the stellar host XUV flux. 

Related t o M D , ISouthworth et all (|2007f ) and 
ISouthworthl (|2008f) show the planet surface gravity, 
g p is also correlated with P. The scatter of the g p -P 
relation is also larger than the observational uncertain- 
ties. XO-5b has a higher g p than the trend based on 
previously known extrasolar planets. It is not clear 
which parameter M p or g D is more fundamenta lly corre- 
lated, if at all, with P. ISouthworth et al.1 (|2007l ) suggests 
the <7 P -P relation is more fundamental than the M p -P 
relation given the direct role g v has on the evaporation 
of hig hly irradiated gas giant planets. iTorres et al.1 
(|2008f ) note that the scatter around their <? P -P trend 
does not correlate with mctallicity. 

The pr eviously known sample of t ransiting planets was 
shown bv lHansen fc Barman! (|2007f) to separate into two 
classes by the Safronov number of the planet. The 
Safronov number, 9 = l/2(V esc /V or b) 2 , where V esc is 
the escape velocity from surface of the planet and V or b 
is the planet orbital velocity, is a measure of the abil- 
ity of a planet to gravita tionally scatter other bodies. 
lHansen fc Barmanl (|2007f ) and refined by ITorres et al.l 
( 20081 ) find transiting planets fall into categories defined 
as Class I, 6 ~ 0.07 ± 0.01 and Class II, 9 ~ 0.04 ± 0.01. 
XO-5b, 6 = 0.1 ± 0.01 does not appear to belong to ei- 
ther class, yet it is not as discrepant as the other outlying 
transiting p lanets, GJ 436b (9 = 0.025) and HAT-P-2b 
(6 = 0.94) ([Torres et al.ll2008j ). There is some overlap, 
but the planets with higher mctallicity stellar hosts tend 
to be the Class II objects with lower 9, whereas Class I 
objects have increasing 9 toward lower metallicities. XO- 
5b does not follow this trend, but has high 9 for higher 
mctallicity. 

A single object such as XO-5b does not invalidate the 
trends amongst the Jupiter-class of transiting planets, 
but illustrates the importance of expanding the number 
of such objects to fully explore the diversity of this popu- 
lation. Recently, the SuperWASP project announced 10 
new transiting planets 11 . We look forward to the detailed 



analysis of these new transiting planets to amplify or di- 
minish the strength of the current trends amongst the 
close-in transiting planets. The end goal of a large pop- 
ulation of transiting planets is to investigate any trends 
that may directly constrain processes that affect planet 
formation, migration, and planet survival. 
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TABLE 1 

XO, E.T., & Perkins Light Curve Data a 



Heliocentric Julian Date 


Light Curve 
[mag] 


Uncertainty 
(l-o-) [mag] 


Filter 


N b 


Observatory 


2452961.13135 


-0.0033 


0.0113 


W 


1 


xo 


2452961.13135 


-0.0095 


0.0113 


W 


1 


xo 


2452961.14526 


0.0286 


0.0113 


W 


1 


xo 


2452961.14526 


0.0048 


0.0113 


VV 


1 


xo 


2452964.11816 


0.0187 


0.0100 


W 


1 


xo 



a The complete version of this table is in the electronic edition of this article. 
The printed edition contains only a sample. 
b Average of N measurements 



TABLE 2 
Stellar Properties 



Parameter 

CSC ID 

RA (J2000.0) 

Dec (J2000.0) 

Galactic Latitude b [deg] 

" Longitude 1 [deg] 

V 

(B-V) 

(V-R c ) 

(V-I c ) 

Vtass 

(V-I c )tass 

J 

(J-H) 
(H-K) 

Spectral Type 
Distance [pc] 
fj, a [mas yr" 1 ] 
/is [mas yr^ 1 ] 
Stellar Mass [M ] 
Stellar Radius [Rq] 

Teff [K] 

[Fe/H] 

\ogg [cm s 2 ] 
logg [cm s~ 2 ] 
v sin i [km s — 1 ] 
[Na/H] 

Stellar p [g cm -3 ] 
Age [Gyr] 

References: 

a) SIMBAD 

b) this work 

c) TASS fDroege et al. 2006) 

d) 2MASS (Skrutskic et al. 2006) 

e) UCAC2 fZacharias eTaT 2004) 

f) SME Spectroscopic Determination (ft I3.lt 

g) Transit Light Curve Determination (§ 13.21 1 



Value of XO-5 R< 



cicrcnce 



02959-00729 
7' l 46 m 51 s .98 
+39°05'40".5 
26.94 
180.63 
12.13±0.03 
0.84±0.05 
0.47±0.05 
0.83±0.05 
12.17±0.06 
0.82±0.08 
10.77±0.02 
0.33±0.03 
0.10±0.03 
G8V 
270±25 
-32.3±2.7 
-24.4±5.6 
1.0±0.03 
1.11±0.09 
5510±44 
0.25±0.03 
4.52±0.06 
4.34±0.07 

1.8±0.5 
0.18±0.03 
1.02±0.2 
8.5±0.8 



a. 
a. 
a 
a 
a 
b 
b 
b 
b 
c 
c 
d 
d 
d 
b 
b 
e 
e 

b.f 

b 

b.f 
b.f 

b,f 

b,g 

b.f 

b.f 

b 

b 



TABLE 3 
XO-5 Radial Velocity Shifts 



Julian Date Radial Velocity Uncertainty 



-245000 


Shift [m s- 1 ] 


[m 


4442.7789 


75.4 


36 


4452.7562 


-147.4 


46 


4456.7451 


-76.6 


28 


4458.9892 


131.2 


20 


4461.7447 


-120.5 


11 


4464.7206 


0.7 


82 


4466.7167 


-24.8 


13 


4476.6851 


110.0 


21 


4478.6743 


-130.2 


11 


4479.6874 


72.8 


24 
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TABLE 4 
Mid- Transit Times 



Heliocentric Julian Date Uncertainty Observatory ID 



-245000U 


(l-a) [day] 




2999.0114 


0.02 


xo 


3354.9557 


0.02 


xo 


3359.1669 


0.02 


xo 


3375.9428 


0.02 


xo 


3417.8079 


0.02 


xo 


3438.7377 


0.02 


xo 


3442.9195 


0.02 


xo 


4129.7088 


0.003 


MF 


4439.6022 


0.003 


EGM 


4464.7268 


0.003 


CF 


4464.7331 


0.003 


CF 


4485.6730 


0.003 


CF 


4485.6685 


0.003 


JG 


4506.6061 


0.003 


MF 


4527.5427 


0.003 


JG 



a Observatory ID is author initials, except XO is XO 
survey data. 



TABLE 5 
The Planet XO-5b 



Parameter Value 

~~P 4.187732±0.00002 d 

t c 2454485.6664±0.0004 (HJD) 

K 145±10 m s- 1 

a 0.0508±0.0005 A.U. 

i 86.8±0.9 deg 

M p 1.15±0.08Mj 

R p 1.15±0.12 Rj 

g p 22±5 m s~ 2 

a/R* 9.81±0.8 

R p /R+ 0.106±0.003 

t 3.13±0.07 hr 

Impact parameter, b 0.55±0.09 

logg p 3.35±0.09 cgs 

p p 1.02±0.3 g cm" 3 

T eq = T eff (R*/2a) 1 / 2 1244±48 K 

Safronov 0.10±0.01 



